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In this article, we propose a mathematical approach that connects an 
innovative spatio-temporal model to the problem of the EG.5 variant of 


COVID-19 in a human population. We demonstrate the existence and 


*Corresponding author 


Received 27 November 2023; revised 19 February 2024; accepted 19 February 2024 
El Mehdi Moumine 


Laboratory of Analysis, Modeling and Simulation, Department of Mathematics and 
Computer Science, Faculty of science Ben M’sik, University Hassan II, Casablanca, Mo- 


rocco. e-mail: moumine.mathsQgmail.com 


Omar Balatif 
Laboratory of Fundamental Mathematics and Their Applications , Department of Math- 
ematics Faculty of Sciences El Jadida, Chouaib Doukkali University, El Jadida, Morocco. 


e-mail: balatif.maths@gmail.com 


Mostafa Rachik 
Laboratory of Analysis, Modeling and Simulation, Department of Mathematics and 
Computer Science, Faculty of science Ben M’sik, University Hassan II, Casablanca, Mo- 


rocco. e-mail: m_rachik@yahoo.fr 


How to cite this article 

Moumine, E.M., Balatif, O. and Rachik, M., Mathematical analysis and forecasting of 
controlled Spatio-temporal dynamics of the EG.5 Virus. Iran. J. Numer. Anal. Optim., 
2024; 14(2): 500-521. https: //doi-org/10.22067 /ijnao.2024.85581.1349 


500 


501 


Mathematical analysis and forecasting of controlled ... 


uniqueness of the global positive solution for our suggested system. The 
implementation and analysis of an applicable optimal control issue are as 
follows. The methods of optimal control theory are applied in this work 
to demonstrate the existence of optimal control, and with necessary op- 
timality conditions, we discover the explicit expression of optimal control 
that minimizes the negative impacts of this infectious disease on countries. 
We provide numerical simulations at the conclusion to demonstrate the 


efficacy of our chosen strategy. 


AMS subject classifications (2020): 49J15, 92B05, 93A30, 93C10 


Keywords: Spatio-temporal Model; EG.5 Virus; Optimal Control; Mathe- 
matical Model. 


1 Introduction 


The COVID-19 pandemic swept into our world suddenly, causing global dis- 
ruption and emphasizing the need for cutting-edge scientific methods to un- 
derstand and stop the spread of the disease. Among these methods, quanti- 
tative analysis and optimal control are essential for understanding the spatio- 
temporal spread of the virus and formulating powerful strategies to fight the 
epidemic. 

Since the discovery of the first cases in December 2019, the epidemic 
has wreaked havoc on the economy, society, and public health worldwide. 
The disease quickly crossed international borders, confounding health agen- 
cies and efficiencies around the world. Scientists and researchers have used 
mathematical techniques to study the dynamics of the spread of COVID-19 
and formulate customized public health policies in this fight against an im- 
placable foe (see [15, 9, 2, 12, 1, 14, 16, 23, 11, 10, 7, 8, 21, 18, 3]). Although 
the pandemic might appear to be a distant recollection for some, the coron- 
avirus persists in its spread, with a surge in new cases attributed to emerging 
mutant strains. One such mutation, known as EG.5, has gained prominence 
in Europe since its initial detection earlier this year (2023). Recently, the 
world health organization (WHO) has classified it as a “variant of interest” 


in response to the escalating global case numbers. 


Tran. J. Numer. Anal. Optim., Vol. 14, No. 2, 2024, pp 500-521 


Moumine, Balatif and Rachik 502 


EG.5 virus represents a subcategory of the Omicron variant of COVID- 
19 and exhibits a close genetic relationship with other globally circulating 
variants. It signifies a modified form of the virus. Its prevalence has surged 
on a global scale, rising from 7.6% of COVID-19 cases in late June to 17.4% 
by the close of July (2023), prompting the WHO to designate it as a matter 
of particular concern. Health authorities worldwide are currently conducting 
investigations into this novel virus. Similar to its predecessors, this fresh 
variant has undergone multiple mutations. Scientists are tirelessly working to 
grasp its attributes, encompassing its transmissibility, severity, and potential 
implications for vaccine efficacy. Presently, the public is strongly encouraged 
to persist in adhering to safety protocols, including mask-wearing, practicing 
social distancing, and receiving vaccinations, in order to mitigate the virus’s 
propagation. 

In addition, it is necessary to model and predict the spread of the virus 
and develop effective control strategies by offering theoretical models that 
take into account the complex interactions between people, the environment, 
and the properties of the virus itself. This mathematical study makes it 
possible to understand the spatio-temporal transmission of viruses. These 
mathematical models, like epidemiological models, provide a framework to 
analyze and predict disease progression in a specific population. Additionally, 
they enable quantification of the impact of various control methods such as 
social isolation, mask use, and vaccination on virus spread. 

However, to maximize the effectiveness of control strategies, an optimal 
control-based strategy must be used. The goal of the latter is to select the 
best preventive measures considering limitations and available resources to 
reduce the spread of disease. This method combines mathematical analysis 
with optimization techniques to pinpoint actions that will have the greatest 
impact on reducing case incidence, mortality and the strain on the healthcare 
system. 

In this article, the optimal management of the spatio-temporal transmis- 
sion of EG.5 virus is reviewed together with current developments in the field 
of mathematical analysis. We look at the spread of the disease as described 
by the epidemiological model. We also move from optimal control strategies 


that we can implement in our spatio-temporal model to create effective pub- 


Tran. J. Numer. Anal. Optim., Vol. 14, No. 2, 2024, pp 500-521 


503 


Mathematical analysis and forecasting of controlled ... 


lic health regulations that take into account the specific characteristics of the 
disease, available resources, and ethical considerations. In order to evaluate 


the outcomes of treatments and to support intelligent decision making. 


2 Background 


However, merely simply simulating how the virus spreads will not be enough 
to stop a pandemic. This is where it is important to have good control. 
The goal of optimal control is to consider the limitations and resources 
at hand and identify proper management methods to reduce virus spread 
[4, 13, 6, 17, 20]. Determining the timing and level of containment mea- 
sures, vaccination campaigns, managing medical resources, or even adapting 
screening and research strategies are some possible solutions. 

The researchers conducted an evaluation of various interventions, employ- 
ing a combination of mathematical analysis and optimal control techniques to 
identify the most effective approaches for reducing viral transmission. In one 
study [15], the authors introduced a model featuring two distinct controls. 
The first control demonstrated an increase in the proportion of individuals 
adopting virus-preventive measures, coupled with a decrease in those not 
doing so. These changes were attributed to media campaigns and direct 
community engagement efforts. The second control illustrated that limiting 
social interactions effectively lowered disease transmission risk, subsequently 
reducing the number of infected individuals. In another study [9], the authors 
aimed to develop a multi-regional model by incorporating four alternative 
strategies as controls. The first control represented awareness campaigns in 
the region “7,” with a focus on disseminating health precautions to protect 
individuals from viral infections. The second control addressed health and 
security measures to restrict movement between regions “j” and “r.” The 
third control involved promoting quarantine measures for exposed individ- 
uals within the same regions, while the fourth aimed to eliminate infected 
animals. 

We will examine a spatio-temporal mathematical model, which is based on 
the one we used in [15]. To gain deeper insights into the virus’s temporal and 


spatial propagation dynamics, our analysis will incorporate a spatial variable. 
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This integration involves introducing a diffusion term utilizing the Laplace 
operator (well explained in [24]). Furthermore, we will conduct an optimal 
control investigation to enhance our understanding. The management of the 
spatio-temporal transmission of EG.5 virus will then be examined, with a 
focus on the role of mathematical analysis and optimal control in evaluating 
the efficacy of current interventions, exploring novel avenues, and offering 
recommendations based on reliable quantitative data. By understanding the 
potential impact of this study, we will be better prepared to meet the future 


challenges linked to this world pandemic. 


3 Model formulation 


The progression of the EG.5 virus in a human population can be explained 


by the SEIHR model that we suggest, which is shown in Figure 1. It is clear 


Figure 1: Diagram depicting the evolution of the EG.5 virus. 


that the EG.5 virus is described in the following temporal dynamic system: 


eo = A — [a E(t) + a21(t)] S(t) — dS (t), 

oo = [a E(t) + aeI(t)] S(t) — (d+ 8) E(t), 

iO) _ ome) — (ara +n) I(t), @) 
uno = nl (t) — (d+ 52 +6) H(t), 

an ~ 6H (t) —dR(t). 


The different compartments of our spatio-temporal model are listed below, 


Table 1, indicating their descriptions (see [15]). It is clear that the temporal 


Tran. J. Numer. Anal. Optim., Vol. 14, No. 2, 2024, pp 500-521 


505 Mathematical analysis and forecasting of controlled ... 


Table 1: The description of the different compartments of our model. 


Compartment Description 


S The number of people likely to be infected by the virus. 
E The number of infected people who are symptomless. 

I The number of infected people showing symptoms 

H The number of hospitalized people 

R The number of the recovered people 


dynamic system (1) cannot adequately describe how the virus spreads over 
the area. We advocate using the Laplacian operator to close this gap. To be 


specific, we add the followings to the deterministic epidemic EG.5 mode: 


=o z) _ dgAS (t,z) +A— (1 —u) [a1 E(t, z) + a2I(t, z)] S(t, z) — dS (t, 2), 

dE u z) _ dp AE (t,z) + (1 — wu) [a1 E(t, z) + aal(t, z)] S(t, z) — (d+ 8) E(t,z), 
ante z) dy; AI (t,z)+ BE (t,z) —(d+61 +n) I (t,z), (2) 
= 2 = dy AH (t,z) +l (t,2) — (d+ 62 +6) H (t,2), 

anti» = drAR(t,z) +H (t,z) —dR(t,z). 


The initial conditions and no-flux boundary conditions are expressed as 
follows: 


S(0,z) = So, E (0,z) = Eo, I(0,z) = Io, H (0,2) = Ho, R(0,z) = Ro, zen. (3) 
Here, 9 represents a bounded domain within R? with smooth boundary 


denoted as 02. One needs to add the equation 


OS OE OL OH OR 
Ov Ov Ov Ov Ou 


0, (t,z) € [0;t,] x an, (4) 


where v represents the outward unit normal vector on the boundary 0Q, and 
ae corresponds to the outward normal derivative on 02. 

Additionally, we have a control function u : [0;t¢] x 2 —> [0;1[, which 
is designed to decrease the count of infected individuals while concurrently 
increasing the count of susceptible individuals. This is achieved through the 


implementation of measures aimed at reducing social interactions, conducting 
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security campaigns to inhibit migration, and promoting health precautions 
to minimize virus transmission. 

Our objective is to optimize both the count of individuals affected by the 
illness and the cost of the treatment program. To achieve this, we formulate 


the problem as follows: Minimize the cost functional represented by 


ty b 
C(S, BL HReu) = ff" at(t,2)dtde + SIult2)iawouy ©) 
Q/0 


subject to the controlled system (2). 
Additionally, we define Ugq as the set of admissible controls, characterized 
by 
Usa = {UE L™(Q),0<u<1 ae. on Q}, (6) 


where Q = [0,t/] x Q, with © representing a bounded region within R? 
featuring a smooth boundary 00. Further details about the parameters of 


our spatio-temporal model can be found in Table 2. 


Table 2: The description of the parameters of our spatio-temporal model (2). 


Symbol Description 


A Recruitment rate 

d Natural mortality rate 

Qy Virus transmission rate between an infected person without symptoms and a person users 
ag Virus transmission rate between an infected person with symptoms and a person users 
B Rate of infected person without symptoms 

Oy Infected viral mortality rate 

b2 Viral mortality rate of an infected person during hospitalization 

n Hospitalized people rate 

@ Recovered rate 

dg Diffusion of susceptible individuals 

dp Diffusion of infected people without symptoms 

dt Diffusion of infected people with symptoms 

dy Diffusion of hospitalized people 

dr Diffusion of recovered people 


4 Existence and uniqueness of global solution 


In this part, we demonstrate that our model (2) has a strong global solution. 
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We consider M(Q) = (L2(Q))”, 
H'(Q) = (ue Lo(Q) : a € L2(Q) and ou € L2(Q) >, 
Ox Oy 
Ou Pu Pui du 
2 _ 1/Q) - 
and H*(Q) = {ueH (Q) : By?’ By?’ Baby’ DyOe 
spaces. 


L2(Q) \ the Hilbert 


Let L? (0,t¢;H?(Q)) be the space of all strongly measurable functions 
v: [0,t¢] —+ H?(Q) such that 


ts 
J let lira < ov. 
0 


Additionally, we define L° (0,t7;H 1(Q)) as the set of all functions v : 
[0,t¢] + H'(Q) that satisfy 


sup (|lv(t, z)lla1@)) < 00. 
te [0,t¢] 


The norm in L® (0,t-; H'(Q)) is given by 
Ill 22 (0, ¢3472(Q)) = inf{me R,: Ilv(t, z)\la2(ay < m}. 


Our spatio-temporal model (2) can be expressed in the following form 


Oy(t, z) 
Ot 


where y= (Yi, Y2sY3> Yas Y5) = Ce 2 Oe ee and - = (fi, fo, fs, fa, fs) is 
defined by 


= Ay(t, z) + f(t, y(t, z)), (7) 


fi =A—(1—) larye + azys] yi — dyn, 
fo = (1 — u) larye + aeys] yi: — (d+ 8) yo, 
fs = By2 — (d+ 41 +) ys, 
fa = ys — (d+ 62 + 8) ya, 
fs = Oya — dys. 

For all k € {1,2,3,4,5} 


6) 
sr = dpAyx + fr(y(t, z)). 


Let A denote the linear operator, mapping from D(A) C M(Q), defined as 
follows: 
Ay = (dgAy1, dpAyp, d;Ays, dn Ays, drAys) 
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y € D(A) = {y = (vis¥esus. ass) € (H?(O))”: 


Oyr — Oy2 — Oyzs — Oya _ OYs 
Ou Ou Ov Ov Ou 


0 on anh. 


Theorem 1. Consider a bounded region 2 in R? with a sufficiently smooth 
boundary. Assuming nonnegativity constraints for y? on (for k € {1, 2,3, 4,5}), 
as well as nonnegative parameters A, d, a 1,2, 3,61,62,7, and 0, and an ad- 
missible control u € Ugqa with an initial condition 

y® € D(A), the system (2) possesses a unique strong nonnegative solution 

y © W'?((0, t¢]; M(Q)), satisfying the following properties: 


Y1; Y2,Y3,Y4, U5 € i (0,¢; H7(Q)) ae Pas (0, t¢; H'(Q)) M L™(Q). 


Moreover, there exists a constant Mo > 0, independent of the control u, such 
that for all t € [0,t,] and all k € {1, 2,3, 4,5} 


Ouk 
Ot 


+ |lYell c2(0,t;,42(@)) + lI¥ellar(ay + lIyelln~¢qy < Mo. (8) 
L?(Q) 


Proof. Since the operator A possesses characteristics of dissipation, self- 
adjointness, and the capacity to generate a Co semigroup of contractions 
on M(Q) [22], it is readily apparent that the function f represented as 
(fi, fe, fs, fa, fs) maintains Lipschitz continuity concerning the variable y 
denoted by (1, ye, 3; Y4; Ys), and this continuity remains uniform with re- 
spect to t within the interval [0,t;]. Consequently, the system possesses a 
unique strong solution y € W!? ([0,t¢]; M(Q)) with 


yx € L? (0, tp; H7(Q)) for all k € {1,2,3,4,5}. (9) 


We prove now that for all k € {1,2,3,4,5}, y, € L°(Q). 


Set c= max { Il fillz=(@)) , Ue nests [ke {1,2,3,4,5}}, and let 
Ug (t, 2) = yx (t, 2) — ct — IIvell n= ca) : 


It is clear that the function U; satisfies the system 
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OU) = d, AU, (t, z) + f(t, y(t, z)) — 6, Be (0, t¢], 


U;,(0, z, 2’) = y2 -— IIvell i ca) ‘ 


This system admits a unique strong solution defined by 


Uelts2) =F(0 (uf llufllamcey) +f PO— 2) (falule)) ~ €) de 


where I(t) is an infinitesimal semigroup related to the operator d,A. 
It comes that U,(t,z) < 0, so yx < ct + IIvelln~ cay: 
Using the similar manner, we can prove that function V;(t, z) = yz(t, z) + 

ct + lIyell n= ca) is nonnegative. So, yz, > —ct — yell L=cay> then |y,(t, z), | < 

ct + [ell Lo cay> which implies that 


yr € L°(Q) V(t,z) € [0,t7] x2, for KE {1,2,3,4,5}. (10) 
Now, we establish that y, € L°° (0, ty; H'(Q)) for all k € {1,2,3,4,5}. Con- 
sider k € {1,2,3,4,5}, and start from the equation 


Oyx (E, Z) 


Fe — dealt, 2) = feltsults2)), (42) € [0,t7] xO 


We have 


[ [ (Aue) (t, z) ~adye(t.2)) ded =f f itive? Sede: 


By employing Green’s formula, we arrive at 


t Oyk 2 t : 
[LE deds-+ di | [ aw) dzds 
o Ja \ ot o Jo 
t OYE t 5 
=2a, f Le — « Aydeds+ [ [ (fr (t, Ye))” dzds 
0 
= ax | ivy} Paz— ak | Vue? a+ f [i (fr ( (t ,Yk)) dzds. 


Since yp € H?(Q) and | yx] ;0(g) ave bounded independently of u, it follows 
that 
ye € L® (0;t¢; H'(Q)), for & € {1,2,3,4, 5}. (11) 


Based on (9), (10), and(11), it is determined that the inequality presented in 
(8) holds. 
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Furthermore, we conclude that the solution (y1, y2, y3, y4; Ys) remains non- 
negative using arguments similar to those employed for the Field-Noyes equa- 


tions in [19]. Consider the set 


T= {(y1, yo, ¥3, ¥4, 5) :0< ye < D fork € {1,2,3,4,5}}, 


and define the convex functions Fy, over II as Fx (yi, y2, 93,44, Y5) = —Ya- 


This leads to the following relationships: 


VU) Flan = VO) Flynn ASO, 

V (2): flyso = V (92): fly-o = — (1 — 4) ooyiys < 0, 
V (Fs): fly=o = V (9s): flys=o = —Py2 S 0, 

V (F4) + flys-0 = V (—y4) + fly=0 = —7Y3 < 9, 

V (F5)- flysso = V (9s) + flys-o = —Oys < 0. 


As shown in [19], the set II is positively invariant. 


5 Existence of an optimal control 


In this part, we aim to establish the presence of an optimal control for mini- 
mizing cost function (5) subject to the reaction diffusion system outlined in 
(2)-(4) with wu € Ugg. The primary outcome we seek to demonstrate in this 


section is as follows. 


Theorem 2. Given the assumptions outlined in Theorem 1, the optimal 


control problem (2)—(5) possesses an optimal solution denoted as (y*,u*). 


Proof. Since u and yx are uniformly bounded in L*(Q) for k € {1, 2,3, 4, 5}, 
the infimum of the cost function is ensured. Let C* = inf,cu,, C(y, u). Con- 
sider a minimizing sequence {un} C Uag such that im C(y", tn) = C, 
where (y?, 93, 93, 2, ye) is the solution of the system (2)—(3) corresponding 


to the control u,. This leads to the following system: 
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a = dgAyf + A— (1 — Un) [ory? + a2y3] yf? — dy?, 

2 = dpAyy + (1 — Un) [arys + ays] yt — (d+ B) yz, 

“us = dyAy8 + By — (d+ +n) yf, (12) 
“rt = dy Ay? + nyR — (d+ 2 + 6) yf, 

Oe = dp Ay + Oy} — dus, 


Oy _ Oys _ Oy3s _ Oy _ Oys 
Ov Ov Ov Ov Ov ct 
Given that H1+(Q)) is compactly embedded in L?(Q), we can conclude 


that y(t, z) is compact in L*(Q) for k € {1,2,3,4, 5}. 


where 


Additionally, we need to establish that {y?(t, z),n > 1} is equicontinuous 
in C ((0, ty], £7()). 


0 n 
The boundedness of st in L?(Q) implies the existence of a positive 


constant A such that 


| [on e2¢:- [on @. 2a: 


<Alt— a. 


Hence, by the Ascoli-Arzela Theorem (see [5]), we can confirm that y7 is 
compact in C ([0,t,],£?(Q)). Consequently, y? > yf uniformly in L?(Q) 
with respect to t. 

Considering the boundedness of Ay? in L?(Q), there exists a subsequence, 
denoted again as Ay?, that converges weakly in L?(Q). Therefore, for all 


distributions y, 


[cau =f wag | ude = f phy. 
Q Q Q Q 
Oyun. OVE 


This implies that Ay? + Ay; weakly in L?(Q). Furthermore, =* > —* 
weakly in L?(Q), and y? + yf weakly in L? (0,t;H?(Q)) and y? > yf 
weakly in L° (0,t¢;H1(Q)). 

Since yy? — yty2 = (yt — yt) v2 + Yi (YF — yz), we deduce that yi? yy > 
yiys in L?(Q). As a result, u, > u* in L?(Q). Furthermore, since Ugg is a 
closed set, it follows that u* € Uaa. 

By taking the limit as n — oo in (12), we can establish that y* is a 


solution to (7) associated with u*. Thus, 
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x Ok “f * b * 2 
Ciutat) = f° aygt,2)dtde + Flu") aren 
0 
ty b 2 
< limint | ays (t, z)dtdz + 5 lu” (¢, 2)llz2(Q) 
0 


ty b 
< tim f ay (t, 2)dtde + 5 |lu"(t,2)Ilz2(@q@) = C*- 
0) 


This confirms that C' reaches its minimum at (y*,u”*). 


6 Necessary optimality conditions 


In this section, we will delve into the optimality conditions for problem (2)— 
(5) with wu € Uga and provide insights into the characterization of the optimal 
control. Consider an optimal pair denoted as (y*,u*), and let u* = u* + eu, 
where € > 0, represents a control function satisfying u € L?(Q) and u € Ua. 
We denote the corresponding solutions of (2) associated with the controls u* 


and u* as y° = (yf, ¥3, 95,44, Y5) and y = (yt, 2, 43,44, Ys), respectively. We 
put 


-d—(1—u*) (a1 y3 + agy3) —(1-4*) ayyt —(1-u*) agut ) ) 
(1—u*) (o1y§ +agu$)  (1—u*)oryt — (d+ 8) (1—u*) ogut 0 0 

G= ty) B —(d +5, +7) 0 0 
ty) 0 n —(d+6g+n) 0 

0 0 0 a) —d 


and 
(ary + O2y3) Ui 
— (ary + 243) Yt 
P= 0 
0 
0 


We have the following theorem. 
Theorem 3. The mapping y : Uaqa > Wt?([0, t¢], M(Q)) for k € {1, 2,3, 4,5} 
is Gateaux differentiable with respect to u*. For any u € Una, yy, (u*) u = Yr, 


consequently, Y = (Yi, Y2, Y3, Y4, Ys) represents the unique solution to the 
problem given by 


ay 
o> AY +GY+uP © subject to Y(0,z) =0. 
E _ yk 
Proof. We define Y¢ = he for k € {1,2,3,4,5}. By subtracting the 


systems corresponding to y; and yz, we derive the following equation: 
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aye € Ee € * € 
Ai = AY*+G*Y* +uP subject to Y*(0,z) =0, for all zE 
with 
—d—(1—u®) (ayy3 + agus) — (1 — u®) ay y§ — (1 — u®) agy§ 0 0 
(1—u®) (ayu3 +oqu§) (1 u®)ayyf — (d+ 8) (1—u®) aguf 0 0 
Ge = 0 B —(d+6, +n) ) i) 
10) 0 n —(d+69+7) 0 
10) 0 0 6 —d 


Let (I'(t),t > 0) be the semi-group generated by A. The solution to this 


system can be expressed as 
t t 
Y*(t,z) = i T(t — «)G*Y* (a, z)dx + i T(t — «)uPdz. 
0 0 


Since the coefficients of the matrix G* are uniformly bounded with respect 
to €, we can apply Grénwall’s inequality to establish that Y; is bounded in 
L?(Q). Consequently, yg > y~ in L?(Q). Taking the limit as « — 0, we 
obtain 


OY 
ot 


=AY+GY-+uP - subjectto Y(0,z)=0, forallze”. 


Applying a similar approach, we conclude that Y; — Y/ as e > 0. 


Consider the adjoint variable q = (q1, 92,93, 44,95), and let G* represent 
the adjoint of the Jacobian matrix G. The dual system associated with our 


problem can be defined as follows: 


2 — Aq— G*q= D* Dw - subject to q(t, z) =9, (13) 
where 
00000 0 
00000 0 
D=}]00100 and y= |a 
00000 
00000 


Lemma 1. Under the conditions of Theorem 1, if (y*,u*) represents an 
optimal pair, then there exists a unique strong solution denoted as q € 
W!*((0, t¢], M@(Q)) for the dual system (13). Specifically, the components 
gz € L? (0,t¢; H?(Q)) NL® (0,t¢; H1(Q)), where k = {1, 2,3, 4, 5}. 
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Proof. This lemma can be established by introducing a change of variable 
t’ = tp — t and employing a similar manner as demonstrated in the proof of 
Theorem 3. 


By employing linkages between the state and adjoint equations, objective 
functional, and standard optimality approaches, we are able to characterize 
the control optimal and get the necessary conditions for the optimal control 


problem. 


Theorem 4. If there exists an optimal control u* and its corresponding so- 
lution y* € W!?((0, t¢]; M/(Q)), then the optimal control u* can be expressed 


as 


u* = min (mms: max (0. (g2 = 1) (ous + anys) Yi )) : (14) 


Proof. Let (y*,u*) be an optimal pair. Let u® = u*+eh € Uga with h € L?(Q) 


and y* be the corresponding state solution. We get 


C' (ut) (h) = lim = (C (u®) - C(w’)) 


= tm? (af f os g yt) dedt+ Ef fw ~ (w’)*) deat) 
= tim (a PLE 8) ae dt 4 At Et 2hu* + ch?) ) ded). 


ww — lim ¥3 ue teh) — Hs 


Since Tim 1 : = Y3, Tim 1y3 = = yz and y5,y3 € 
eS 
L™(Q), it isitewe that C is Gateaux differentiable its respect to u*, and 


we obtain 
tf tf 
= | | aY3dzdt + vf | hu*dzdt 
0 Q 0 Q 
ty 1 
0 0 


If we put h = v — u*, then 


ty 1 
C" (u*) (vu —u*) =| (De, DYyat-+ | (bu*,v — u*) pacqy dt. 
) 0 


Since 
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ty ts 
| (De, Dy at = | (D* Dy, Y) dt 
0 0 
t / Oq ‘ 
-{ (St — 4g— Gray ) at 
‘6 / ay 
= [lav —cr)a 
ts 
= fo aP(v-w)yae 
0 


and Uga is convex. Therefore, C’ (u) (v — u) > 0 for every v € Uaa, which is 
equivalent to 
- (P*q + bu*,v — u*) pag) dt 2 0 for every v € Ua. 


(q2 — M1) (a1y3 + a2y3) yi 
b 


So, bu* = —P*q, thus u* = 
we have that (14) holds. 


. Since u* € Uaa, 


7 Numerical simulation and presentation of the results 


To achieve the interest of optimal control and evaluate the precision of our 
spatio-temporal model with and without control, we carried out a numerical 
simulation of our optimal control problem in the context of modeling the 
spread of disease. We adopt a rigorous methodological approach. Our op- 
timality system is formulated by state equations with initial and boundary 
conditions (2)—(4), adjoint equations with transversality conditions (13), and 
a characterization of optimal control (14). To solve this system, we use the 
forward-backward scanning method in an iterative process. The temporal 
discretization of the state equations is carried out using the explicit Euler 
method, while the adjoint equations are solved retroactively in time. The 
space considered, representing a city Q = 50Km?, is initially (t = 1 day) 
populated randomly in the city, and it is assumed that the epidemic began 
to spread randomly in different areas of the city. The results will be presented 


as follows. 


The values of the parameters connected to our spatio-temporal model 


with and without control are presented in Table 3. 
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Table 3: Parameters for the model without and with control. 


Parameters Values 
(A, d, a1, 02, 8,61, 62,7,0) (0.06,0.125,0.65,0.65,0.5,0.0001,0.001,0.05,0.02) 
(ds, dz,dz,dxy,dp) (0.2,0.2,0.2,0.2,0.2) 


Exposed -t=10 


10 20 30 40 50 


Figure 2: The behavior of compartments FE and I without control 


Before the enforcement of control measures, a pervasive presence of 
asymptomatic infections was observed in multiple locations throughout the 
region, suggesting a widespread dissemination of the epidemic within the city. 
Figure 4 illustrates a persistent high percentage of asymptomatic cases, con- 
sistently standing at 11% of the total population over time. This phenomenon 
played a role in the formation of a substantial cluster of symptomatic indi- 
viduals across the city, comprising 31% of the total population. The visual 
representation underscored the extensive transmission of the infection in var- 


ious areas of the city. 
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Figure 3: The behavior of compartments F and J with control 


After the enforcement of stringent control measures, a noticeable trans- 
formation has taken place in the landscape of the epidemic. These measures, 
which strategically aimed to curtail the widespread prevalence of asymp- 
tomatic infections in various locations across the region, including crucial 
components such as reducing social interactions, conducting safety campaigns 
to discourage migration, and promoting health precautions to minimize virus 
transmission. The proactive interventions have resulted in a significant de- 
crease in the overall percentage of asymptomatic cases, plummeting from the 
initial 11% to a mere 0.2%. Consequently, there has been a substantial reduc- 
tion in the number of symptomatic individuals, with the affected population 
now comprising only 0.3% of the total. This positive development indicates 
the effectiveness of the implemented controls in containing the spread of the 
disease and alleviating its impact on the population. The newly established 
situation reflects a more controlled and manageable scenario, underscoring 
the success of the measures in restricting the transmission of the infection 


within the city. 


Percentage of | (Total Population = 1250) 
100 I 


Percentage of S (Total Population = 1250) Percentage of E (Total Population = 1250) 


= = “with Contr 


Percentage 


Figure 4: The percentage distribution of compartments S, E, and J relative to the total 


population 
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8 Conclusion 


The spatio-temporal model used in this work captures population size, move- 
ment models, and social relationships that are important in the evolution of 
EG.5 virus transmission, and as our study suggests, we now understand the 
role and potential of controls using intervention strategies, such as advocacy 
campaigns, creating social barriers, and focused testing on, to prevent the 
spread of the virus while reducing associated costs and societal disorders. 

In conclusion, our analysis of this spatio-temporal epidemiological model 
contributes to well-known decision analyses and important mechanisms of 
epidemic development. This approach has the potential to combine rigorous 
mathematical analysis with real data to determine effective targeted ther- 
apies, ultimately supporting global efforts to combat the virus and protect 
humanity wall. Finally, to be more accurate and effective models, we can im- 
prove these models by combining other factors such as age, socio-economic 
characteristics and behavioral models. Additionally, efforts should be made 
to improve the data collection system and to develop reports that will pro- 
vide modern, reliable data for accurate sampling, so that we can develop a 


reduction strategy epidemic and protect human health. 
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